Nonlinear modal interactions in clamped- clamped mechanical resonators 
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A theoretical and experimental investigation is presented on the intermodal coupling between the 
flexural vibration modes of a single clamped-clamped beam. Nonlinear coupling allows an arbitrary 
flexural mode to be used as a self-detector for the amplitude of another mode, presenting a method to 
measure the energy stored in a specific resonance mode. Experimentally observed complex nonlinear 
dynamics of the coupled modes are quantitatively captured by a model which couples the modes 
via the beam extension; the same mechanism is responsible for the well-known Duffing nonlinearity 
in clamped-clamped beams. 
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An important topic in nanomechanics is the motion 
detection of mechanical resonators. Several schemes 
have been proposed to attain sensitivities near the quan- 
tum limit of mechanical motion pQ, whereas application- 
driven research is focussed on on-chip detection [2] and 
readout of resonator arrays [3]. Central in any detec- 
tion scheme is the coupling of a mechanical resonator to 
another system, which transduces the motion into a mea- 
surable quantity. Examples of sensitive detectors include 
a single-electron transistor [4], a microwave cavity [5], 
or an optical interferometer [6]. A second mechanical 
resonator can also be used to detect the motion of the 
resonator [8]. Such a system of coupled resonators 
has been proposed as a quantum nondemolition detection 
scheme, in which one resonator is in a quantum state [9]. 
Coupling between different mechanical resonators is of- 
ten present in large-scale integrated arrays due to elec- 
trostatic [7 and mechanical interaction [8 . Coupling be- 
tween individual resonators can also lead to complex be- 
havior, which is theoretically well-documented [ID] . 

In this Letter, we study the coupling between vibra- 
tional modes in a single beam resonator. We demonstrate 
that flexural modes are coupled by the displacement- 
induced tension in the beam. Using this coupling, the 
displacement of any mode can be detected by measuring 
the response of another mode, making otherwise unde- 
tectable modes visible. We present a general theoreti- 
cal framework based on the Euler-Bernoulli equation ex- 
tended with displacement-induced tension. The model 
quantitatively describes the complex dynamic behavior 
observed in the regime where two modes are simultane- 
ously driven nonlinear. The coupling mechanism plays 
an prominent role in the dynamics of carbon nanotube 
resonators and resonators under high tension, and should 
be taken into account when describing such systems ac- 
curately. 

Experiments are performed on a single-crystalline 
silicon beam with dimensions L x w x h = 1000 x 
35 x 6 jam 3 fabricated by patterning a silicon-on-insulator 
wafer and subsequent wet etching. The resonator is 
placed in a magnetic field of B = 2.1 T and a magne- 



tomotive technique [3l [TTj is used to detect the mechan- 
ical motion of the beam at room temperature and at- 
mospheric pressure (see Figure [T^i). The beam is driven 
at multiple frequencies by sending alternating currents 
through a conductive aluminum path, evaporated on top 
of the resonator. The motion of the beam in the magnetic 
field generates an electromotive force, which is balanced 
using a Wheatstone bridge, amplified, and then digitized 
(see Figure |TJl) . The frequency response (amplitude and 
phase) of the resonator at the two drive frequencies is 
calculated using a digital signal processor. 

Measurements are conducted on the first and third 
flexural mode of the beam. In Fig. [T]o the frequency re- 
sponse of the third mode is shown for three different drive 
amplitudes of the first mode. The resonance frequency of 
the third mode increases when the amplitude of the first 
mode | Ai| becomes larger. At the same time the am- 
plitude decreases slightly. The resonance frequency /r,3, 
obtained by fitting a damped driven harmonic oscillator 
response, is plotted as a function of the drive amplitude 
of the first mode in Fig. For small amplitudes 
a quadratic dependence of /r } 3 on is found (dashed 
line in Fig. [1^). To give a qualitative picture, which ex- 
plains the observed coupling in Figure ^p. When driving 
the third mode on resonance, the tension in the beam 
increases as the amplitude of the first mode increases. 
This results in a higher frequency and a lower amplitude 
of the third mode. 

As mentioned above, using the magnetomotive mea- 
surement scheme it is not possible to detect the even reso- 
nance modes. However, by exploiting the coupling to de- 
tectable odd modes, their motion is observable. We have 
detected the motion of the second mode by measuring its 
influence on the first mode. To drive the second mode, 
the sample (different device, h = 10 jam) is mounted on a 
piezo actuator and excited at frequencies around the sec- 
ond resonance mode. Off-resonance, no frequency shift of 
the first mode is observed. When driven at the resonance 
frequency of the second mode, the resonance frequency 
of the first mode shifts to a higher value (Fig. The 
response of the second mode is obtained by measuring 
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FIG. 1: (a) Setup, with a colored scanning electron micro- 
graph of the resonator beam. DSP is the digital signal pro- 
cessor, R is used to balance the bridge and set to 22 Q and A is 
the amplifier, (b) Frequency responses of the third mode (am- 
plitude | A3 1) for different drive amplitudes of the first mode 
I Ai| on resonance. The drive current of mode 3 is I3 = 1.5 
mA. Inset are the beam shapes of the first (solid) and third 
mode (dashed) for increasing amplitude of the first mode. 
For clarity, the amplitude of the third mode has been lowered 
more than in the actual situation, (c) Resonance frequency of 
mode 3 (/r,3) as a function of drive amplitude of mode 1 with 
measured data (black squares) and result from the model (red 
line). For small drive currents a quadratic dependence is ob- 
served (dashed blue line). Inset: schematic of the resonator 
in a magnetic field showing the beam shape of the first (blue) 
and third (red) mode, (d) On a different device, the resonance 
frequency of the first mode is used as a detector of the second 
mode. Inset: beam shape of the first (blue) and second (red) 
mode. The error bars are within in data markers. 



the shift in /r 5 i, whereas a magnetomotive measurement 
around the same frequency shows no signature of the sec- 
ond mode. For weak piezo driving the frequency shift is 
proportional to | ^4_2 1 2 and a squared damped driven har- 
monic oscillator function fits the data well with fi = 287 
kHz and Q = 250. 

In order to quantify the coupling between the flexu- 
ral modes of the beam, an analytical model is developed. 
First, the equations are derived for the general situation 
with modes coupled; then we focus on the experimen- 
tal situation, where only two modes are considered. The 
Euler-Bernoulli equation including tension T }12HT4] is 
used as starting point. To simplify the notation, the dis- 
placement u is scaled with the beam thickness /i, and the 
coordinate x with the beam length L. A displacement of 
the beam causes an elongation and increases the tension. 
The dimensionless tension T = L 2 T/D (D is the bending 
rigidity) is given by: 
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To is the residual tension in the beam and r = h 2 A/I y , 
with I y the second moment of inertia and A the cross- 



section. For a rectangular beam r equals 12 [15 . The 
displacement u and dimensionless force F = L A F/Dh 
can be split into a dc and an ac part i.e., u = ua c + u ac 
and F = F^ c + F ac . This yields a well-known equation 
for the static displacement [T2HT4] and an equation for 
the ac motion: 

^ac+^ac+^M-^-Tde)^^-^-?-^)^ = T ac . 

(2) 

Here, T& c is the residual tension plus the tension from 
the dc displacement and T is the time-averaged tension, 
which also contains terms proportional to u^ c . For small 
ac displacements T « T dc . T ac contains all terms that 
are linear in u ac . The operator C[u) is defined as [T4] : 



C[u) = u"" -T dc it" -T ac [u}u[ 



dc- 



(3) 



The first three terms on the left side of Eq. [2] determine 
the linear response of the system. The nonlinearity is 
introduced with u 2 c and u\ c1 which occur in the last two 
terms. 

The resonance frequency for infinitesimal small am- 
plitudes cJq^ and the corresponding shape £i(x) of mode 
i are the eigenvalues and orthonormal eigenfunctions of 
C respectively. The ac displacement is expanded in 
terms of the mode shapes as u ac = Y^=i u i{^)^i{ x )- 
The dc displacement is defined as i^c = ^o^o(^) and 
Jq 1 ^ , i (x)^j(x)dx is denoted as Iij. The value of the inte- 
gral Iij depends only on the shapes of mode i and j and 
and can be calculated numerically. Using this notation 
the tension T is: 



T(t) 



v*{t)uj{t)Iij> (4) 
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so that T dc = T + t^5%o and T ac = \tu q Y^=i u n I n o- 
For nonzero amplitudes u n the tension increases, tuning 
the resonance frequencies uj^i away from ujq^. The effect 
of the resonance frequency on its own motion results in a 
Duffing equation (i.e., the resonance frequency increases 
with its own amplitude). Based on the same concept, the 
tuning of the resonance frequency due to the motion of 
other modes can be envisaged and this coupling is the 
central theme of this work. 

The displacement of a beam, which is driven at fre- 
quencies uoi is written as: 



u(x, t) 
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ai\£i(x)cos(u)it + Zcii), 



(5) 



where the = Ai/h are the complex amplitudes of the 
mode at uji. Substituting Eq. [5] into Eq. [4] gives an ex- 
pression for the total tension. The two tension terms in 
the ac equation are given by: 
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T-T*c = ^(^|^|% + y^^ + cx.)+(T-T-T ac ) 

i>0 

T-T-T ac = - {^-a\j 2 ^Hii + a^ (6) 

i#i>0 



Here, c.c. stands for the complex conjugate of the term 
before c.c. between the round brackets. The expressions 
for the tension, Eq. [6j are substituted in the equation 



where F ac = L^T^jDh is the dimensionless ac force. 

So far, the analysis is valid for any flexural res- 
onator. We now focus on the experimental situation 
where, unlike in buckled beams [13] and string-like res- 
onators [16], the residual tension does not play a signif- 
icant role, To = 0. Moreover, the static displacement 
uq = 0, which may not be the case for carbon nanotube 
resonators [17], where a gate voltage induces a static dis- 
placement. We finally assume that the resonances are 
resolvable: \u R ^ - u RJ \ > u Rii /Qi + urj/Qj for i ^ j. 
The homogeneous force per unit length on the beam is 
/ac = BI, where / is the current through the resonator. 

To compare the model with the data in Fig. [TJ^, we 
first extract the experimental values of the parameters. 
The resonance frequencies for the first and third mode 
are fi = 48.2 kHz and f 3 = 273.4 kHz, close to the 
predicted values of 46.1 and 249 kHz for a beam-like res- 
onator. Their Q-factors are Qi = 41 and Qs = 172. The 
values of Iij determine the coupling strength: In = 12.3, 
I 33 = 98.9 and /13 = I31 = —9.7. The average displace- 
ment of the modes per unit deflection, J Q ^i(x)dx^ are 
0.83 and 0.36 for mode 1 and 3 respectively. The model 
is solved numerically by calculating the amplitudes of 
the two modes self-consistently for the experimental con- 
ditions and without any free parameters. The calculated 
resonance frequency of the third mode as a function of 
the amplitude of the first mode is shown in Fig.[T]3. Excel- 
lent agreement is found between the observed frequency 
shift and the prediction by the model. The proposed 
model which couples modes though tension induced by 
the vibrations is therefore a good description of the ex- 
periment. For large amplitudes, the resonance frequency 
scales with |Ai| 2 / 3 , indicating that the beam is in the 
strong bending regime [12]. For small the tuning is 



I 

of motion, Eq. [2] The time-averaged equation for the 
amplitude of mode i driven at frequency u^, is then given 
by: 



(7) 

I 

quadratic and the third mode can be used to detect the 
amplitude of the first mode with a sensitivity of 0.18 Hz 
nm -2 , which is determined from the quadratic curve in 
Fig. [5. 

To further test the consistency of our model, we 
study the complex dynamics of the coupled modes. When 
driving both modes nonlinear, interesting features are ob- 
served. In Fig. [2^i the amplitude of the nonlinear first 
mode is plotted versus the driving frequencies fi and 
f%. Simultaneously the amplitude of the third mode is 
recorded (Fig. [2]o) . The two modes interact with each 
other as the nonlinear line shape of one mode is reflected 
in the response of the other mode. Also a frequency re- 
sponse with two peaks, which is clearly different from a 
Duffing line shape, is observed as illustrated more clearly 
in Fig. [2^ and f. The two peaks arise from the bistable 
first mode, where two values for the amplitude are pos- 
sible. These two amplitudes correspond to two values 
of the tension, which leads to two resonance frequencies 
of the third mode and to two peaks in its frequency re- 
sponse. The simulation with the parameters as stated 
above, reproduces all observed features in the amplitude 
of both modes (Fig.^-d). This indicates that the model 
captures the coupling mechanism in detail [15] . 

An example of how the coupling between the modes 
can be used in practice, is the increase in dynamic range 
of the mechanical resonator. In small-scale resonators 
the dynamic range is limited by the nonlinear response 
at strong driving amplitudes, which is disadvantageous 
for many applications [19] [20] . Our analysis and exper- 
iments show that there is a way to extend the dynamic 
range of one mode by driving another mode on resonance 
at high amplitudes. Fig. [3^i,c shows the frequency re- 
sponse of the third mode when the frequency of the first 
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FIG. 2: Frequency- frequency response of the simultaneously 
driven and detected first and third mode. The drive frequency 
of the third mode is swept, while driving the first mode at a 
fixed frequency. After each sweep the frequency of the first 
mode is increased. The amplitudes \Ai\ (a) and |^4_3 1 (b) are 
recorded. The driving currents are ii = 0.8 mA and ^3 = 7.0 
mA. Red indicates a high amplitude and blue corresponds to 
a low-amplitude response. Simulations are shown in (c) and 
(d). The frequency response is plotted for the amplitude of 
mode 1 (e) and mode 3 (f) at fi = 48 kHz, also the back 
sweep is shown for the two cases. 



mode is swept across its nonlinear resonance. Away from 
/r ? i the third mode shows a Dumng-like response as il- 
lustrated in the top and bottom panel of Fig.[3]3,d. How- 
ever, when driving the first mode on resonance (middle 
panel) the third mode displays a hysteresis- free response. 
Effectively, the nonlinearity constant in the Duffing equa- 
tion is decreased, which can be understood as follows: 
when the third mode enters its resonance, its amplitude 
increases and the increased tension tunes the resonance 
frequency of the first mode up. The amplitude of the first 
mode then drops, reducing the tension and lowering the 
resonance frequency of the third mode. This feedback 
mechanism reduces the cubic stiffness of the third mode 
and makes the third mode linear, thereby increasing the 
dynamic range. 

The presented model applies to any clamped- 
clamped geometry, ranging from suspended bridges to 
carbon nanotubes. For nanomechanical devices with high 
aspect ratios, the flexural rigidity can be neglected and 
the restoring force originates from the axial rigidity. The 
energy of the resonator is stored in the beam elongation 
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FIG. 3: (a) Frequency- frequency plot for the forward sweep 
of the third mode. Red indicates a high amplitude and blue 
corresponds to a low amplitude. On resonance of the first 
mode, the third mode is linear and off-resonance it is nonlin- 
ear, which demonstrates the tuning of the nonlinearity. (b) 
Forward (black lines) and back (red lines) sweeps are shown 
below resonance, on resonance and above resonance, with 
fi = 45, 52 and 55 kHz respectively. Driving currents are 
ii = 2 mA and 7 3 = 8 mA. (c) Result of the simulation with 
parameters extracted from the experiment, (d) Simulated 
traces for the situation in (b). 



which couples directly to the ac tension. This increases 
the coupling between the resonance modes and makes 
the detection mechanism well-suited for nanowires and 
nanotubes. Furthermore in nanomechanical devices, sig- 
nificant residual tension (To) may be present in the beam. 
Only the numerical values 1^ change in that case. For a 
resonator with high To, the off-diagonal elements Iij = 
and for mode 1 and 3, In is 9.9 and I33 is 88.8. Thus, 
for a string-like resonator the coupling term in Equa- 
tion [7] remains within the same order of magnitude as in 
the beam resonator used in the experiment. To quan- 
tify the effect of modal interactions in a nanomechanical 
resonator with tension, we consider the suspended car- 
bon nanotube with a high-quality factor from Ref. [2T] . 
Taken the parameters listed in that paper, we calculate 
the sensitivity of the third mode to the amplitude of the 
first mode. We find a value of 1 MHz nm -2 , which is 
more than six orders of magnitude larger than the value 
found from the dashed line in Fig. 

In conclusion, we have measured the coupling be- 
tween flexural modes of a clamped-clamped beam res- 
onator by simultaneously driving the beam at multiple 
frequencies. We observe nonlinear interaction between 
the modes in the linear regime, and complex dynamics 
at large driving amplitudes. When describing the motion 
of a mechanical resonator, it is necessary to include this 
interaction, since this mechanism divides the available 
energy over the modes, and plays a role in the energy 
dissipation in the resonator. A theoretical model is de- 
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veloped, which couples arbitrary flexural modes via the 
tension. The model is in excellent agreement with the 
measurements and quantitatively captures the observed 
complex dynamics. The nonlinear coupling can be used 
to detect resonance modes that would otherwise be in- 
accessible by the experiment, to tune the nonlinearity 
constant, and to increase the dynamic range of micro- 
and nanomechanical resonators. 
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